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Abstract 

There is renewed interest in the question of whether 
the Navier-Stokes equations (NSE) , one of the fun- 
damental models of classical physics and widely 
used in engineering applications, are actually self- 
consistent. After recalling the essential physical as- 
sumptions inherent in the NSE, the notion of weak 
solutions, possible implications for the energy con- 
servation law, as well as existence and uniqueness 
in the incompressible case are discussed. Emphasis 
will be placed on the possibility of finite time sin- 
gularities and their consequences for length scales 
which should be consistent with the continuum hy- 
pothesis. 



1 Introduction 

As computational fluid dynamics makes progress to- 
wards the simulation of realistic three-dimensional 
flows, the validity of the Navier-Stokes equations 
(NSE) can be tested in a more and more refined 
way. To put it from an applied point of view: Be- 
fore experiments in wind tunnels are substituted by 
computer simulations, one should make sure that 
the underlying theory is at least self-consistent. As 
a matter of fact, after the classical mathematical 
work by Leray Hopf , Ladyzhenskaya , Ser- 
rin [^, Temam to refer to important contribu- 
tions in the field, there is renewed interest in the 
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fundamentals of the NSE, see for instance the mono- 
graph of Doering and Gibbon [Q, or a series of pa- 
pers by Lions Q and references therein. 

This contribution focuses on the question of self- 
consistency which arises, when one of the assump- 
tions inherent to the NSE, namely the continuum 
hypothesis, is confronted with the length scales 
emerging from solutions of the deterministic NSE. 
First, the NSE will be briefly derived from physical 
principles with due attention paid to the continuum 
hypothesis. After recalling the notion of weak solu- 
tions, the state of the art of mathematical existence 
and uniqueness proofs will be indicated. The impli- 
cation of weak solutions upon energy conservation 
will be discussed. The possibility of finite time sin- 
gularities will be related to length scales and thus 
to the problem of self-consistency. 

2 Derivation of the NSE 

The NSE are based on the conservation of mass and 
on Newton's second law. In addition, the more spe- 
cial assumption of a so-called Newtonian fluid is 
adopted, which is justified in a great many cases 
of hydrodynamic fiows. To formulate the conserva- 
tion laws it is customary to pick out a connected 
cluster of molecules contained in volume Vt which 
is deformed in time and translated according to the 
local velocity v(x, t) of the fiow. Time derivatives 
of corresponding magnitudes are conveniently eval- 
uated by means of the Reynolds transport theorem 

^ dVfix,t) 



dV 



dV 



dt jVt 
df 



dt 



+ div(/v) 



(1) 



1 



where / is a scalar function. If p(x, t) denotes the 
mass density, then conservation of mass, namely 



dt Ivt ^ 



gives rise to the continuity equation 
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which if p = po is constant, leads to the incompress- 
ibility condition 



div(v) = 0. 



(4) 



Newton's second law implies that any change in mo- 
mentum is caused by external forces which in con- 
tinuum physics are described by the volume force 
density f (e.g. gravity) and by a tensorial force 11. 
This tensor reflects the influence of the adjacent 
fluid on a given fluid particle. The momentum bal- 
ance reads 
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where n dS is the oriented surface element of the 
volume Vj. It is convenient to separate in 11 an 
isotropic part, the pressure p, which is present also 
in the hydrostatic case, from the so-called stress 
tensor T 
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-pSik + Tik, i,k = 1,2,3. (6) 



The Newtonian fluid assumption now amounts to 
the following linear relations between T and the 
strain (rate) tensor S: 
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The 4th rank tensor C is constant and describes the 
effect of viscosity. In the isotropic case, C is of the 
form 



Cikmn — ^^ik^mn "I" l^i^im^kn "J" ^in^kn 



(9) 



where ji and v are macroscopic viscosity parame- 
ters. In the incompressible case, u drops out, and 
after making use of mass conservation we can write 
down the momentum balance as follows 
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Here Vt is an arbitrary local space volume. To be 
sure of the existence of the above integral, one may 
adopt the sufficient conditions that the following 
fields are locally square integrable 
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This can be easily seen with the aid of the Schwarz 
inequality. For instance, if Xj is a cartesian unit 
vector, then we can write 
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As will be discussed later on, there may arise dif- 
ficulties with the conservation laws when certain 
weak conditions on the velocity field v are adopted 
as is customary in the frame of functional analysis. 
From eq.(lO), the following standard NSE in the 
form of partial differential equations are inferred 
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= -grad(p) + /xAv + /9of (13) 



where p is determined through the incompressibility 
condition div(v) = 0. 



3 Continuum assumptions 
length scales 



and 



The NSE describe macroscopic physical quantities 
which constitute mean values with respect to the 
underlying atomic degrees of freedom. The den- 
sity p(x, t) at the space point x, for instance, has 
to be understood as an average over some volume 
AV centered at x. If AV is chosen too small, a 
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single measurement of p may largely deviate from 
its mean value due to molecular fluctuations. An 
estimate for a physically reasonable lower bound of 
AV^ can be deduced from the mean thermal den- 
sity fluctuation A/o as given in standard textbooks 
of thermodynamics H 



(14) 



where k is the Boltzmann constant, T the absolute 
temperature and k the compressibility. If we require 
the relative fluctuation ^pj p to be smaller than, 
say 10^3, at T = 300° Kelvin, then we find that 
the diameter d of the volume Ay should be d > 

3 • 10~''m for air, or d > 10~^m for water. 

As an implication, if solutions of the deterministic 
NSE turn out to vary on a space scale much smaller 
than the above lower bounds, then we are outside 
of the validity domain of these equations. Here is 
the point where the self-consistency problem arises. 

In the turbulent regime, length scales decrease 
with increasing Reynolds number R. As is listed 
e.g. in Q, the Kolmogorov length bx below which 
eddies are destroyed by dissipation, is given by 
bK = LjW'l^ where L is a typical external length, 
e.g. the diameter of the containment. As another 
example the thickness (5^ of a turbulent boundary 
layer scales as bs ~ L/(RlogR). If L = 1 cm, then 
bx and bB reach the continuum limit at i? ~ 10^. 

4 Weak solutions and energy bal- 
ance 

Since Leray's pioneering work one has been 
looking for generalized solutions v(x, t) of the in- 
compressible NSE in the space time domain ■= 
ri X [0, r] with the following properties: 



v(x,t = 0) = q;(x) 
v|af7 = 0, 
div(v) = 0. 
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The above equations correspond to the initial con- 
dition, no-slip boundary condition and incompress- 
ibility, respectively. To establish weak solutions. 



test vector fields $ € S are introduced with the 
following properties 
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where D(i}) is the Schwartz space (C°° and com- 
pact support in 17). Now v(x, t) is called a weak 
solution if it is locally square integrable and if the 
following projections of the NSE and the continu- 
ity equation hold for every <I> G D and for every 
scalar function (p with compact support in 17, 
respectively j0] 
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where z/q = p/po denotes the kinematic viscosity 
and summation convention is adopted. The pres- 
sure term dropped out in (19) due to the solenoidal 
property of <I>. A typical theorem reads |Q: 

Theorem: A unique weak solution exists, at least 
in the time interval t G [0, ri] with ri < r, pro- 
vided the initial velocity field Q:(x) G 1^2^^ and the 
external force density f obeys the condition 
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where denotes the Sobolev space with the sec- 
ond space derivatives being square integrable. 



As should be noticed, even if the condition ( |2l| ) 
on the external field f holds for arbitrarily large r, 
uniqueness can be guaranteed by the above theorem 
only for the smaller time interval t G [0,ri]. While 
this is typical in space dimension three, one has 
Ti = r in the case of two-dimensional flows. 

Which price do we have to pay for accepting weak 
solutions? To discuss a possible implication for en- 
ergy conservation, we recall the notion of weak and 
strong convergence of a sequence of real functions 
a^^\a^'^\ ..a^^\ ... This sequence is called to con- 
verge weakly against the function a*, if for any 
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square integrable function g 

dya(^)a(^) <oo and 

lim / dVa^'^^g= [ dVa*g. (22) 

It converges strongly, if 

lim / (iya(^)a(^)= / dVa*a*. (23) 

In the case of weak convergence, we have the iden- 
tity § 

lim { / dV (a(^) - a*f - [ a^^^a^^) + 



for the weak solution v* together with the test field 
$ = v(^) E D and obtain 
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/ dVa*a*} = (24) 
Jo. 



Apart from partial integrations, we made use of the 
incompressibility condition (20) which implies the 
relation 



which is true because / dV a^^^a* converges 
(weakly) to / dV a*a* and the two non-converging 
terms J^dV cancel each other identically. 

As a consequence one has in particular 

liminf / dVa^^^a^^^ > f dVaa*. (25) 

Here, the equality sign is guaranteed only in the 
case of strong convergence where simultaneously 
liminf = lim sup. 

To derive the energy balance for a sequence of 
approximations v*^^^ which converge weakly against 
a solution v* of (19), we use basis functions ^^'^^ € 
D with the properties (18) as 

N 

vW(x,t):=5]c(,)(t)cDM(x), C(,)GR. (26) 

It is convenient to introduce the following abbrevi- 
ation for the kinetic energy at time t 

E(^Ht):=l[ dy.f)(x,t).f)(x,t). (27) 
2 Jn 

E* [t) denotes the energy corresponding to the weak 
solution 

V* = v(^) + r(^) (28) 

where r^-^^ is the remainder to the approximate field 
v(^). We now insert into (19) the above expression 
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The above surface integral vanishes because E 
D. It should be noticed that eq.(29) holds true for 
any cutoff N; it follows strictly from the definition 
(19) of a weak solution; in particular, no approxi- 
mate projection scheme was adopted as is common 
in Galerkin representations. 

In the case of strong solutions with v* € W2 
in the space time domain fi,-, one can show that 

R{N) 

— !■ in the limit ^ 00 so that we would 
have the physically plausible energy balance 

JO Jn OXi dxi 

= E*{Q)+ r dt f dVvlfk, (32) 
Jo Jn 

or in words: the kinetic energy at time r plus the 
energy dissipated up to r equals the initial kinetic 
energy plus the work done by the volume force f up 
to time T. 

However, if v* is a weak solution, then we have 
only the property of boundedness of the integrals 
in (29) and (30), except for the /^-integrals and the 
initial energy EW(0) which converges under the 
assumptions specified in Theorem A. Making use of 
the inequality (25) we can write 
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E*{t) + U0 dt dV^-^ + L, (33) 

Jo JQ OXi dXi 

where L,, > 0. With R^, denoting the hmes inferior 
of R^^\ the energy balance (29) reads in the same 
limit 

E'ir) + .„ rat fav'^'^ 

Jo Jn OXi dxi 

= E*{0)+ r dt f dVvlfk + R^-L^. (34) 
Jo Jn 

Thus, in the case of weak solutions there may be 
unphysical sources or sinks (depending on the sign 
of — L^) of kinetic energy due to the presence 
of singularities. The latter are connected with the 
space gradients of the velocity field, since 
t € (0, r) can be shown to converge under rather 
general assumptions, see also If R^ — < 0, 
then the kinetic energy E*{t) is smaller than phys- 
ically expected; this is known as Leray inequality, 
see e.g. p. 104 of 

5 Uniqueness and finite time sin- 
gularities 

As already mentioned, one gets square integrable 
solutions V under rather general assumptions on the 
external data. The main basic problem of the NSE 
is related to uniqueness which so far is tied to the 
existence of the following time integral, for a recent 
discussion see [01, 



/(r) := / dt \\Dv\ 



with the supremum norm 
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(36) 



The origin of this integral will be indicated in the 
Appendix. Up to now, in three dimensions the ex- 
istence of I(t) has been corroborated only for finite 
time intervals r. If /(r) exists for arbitrarily large 
r, then both uniqueness and existence of weak solu- 
tions can be established for arbitrarily large times 
under quite general conditions. 



If /(r) exists only up to some time r*, then 
IlL'vlloo is singular at t = r* in a way, that there 
is at least one space point xq € il, where one of 
the components dvi/dx^ diverges, for instance as 
follows 



5fi(x, t) 



dxt 



a 



t<T*, 7>1. 



(37) 

Since for t near r* the behaviour (|37| ) implies 
changes of the velocity field over arbitrarily small 
length scales, it is in conflict with the continuum 
assumption. The length scales are then small com- 
pared to the diameter of the volume AV of a 
fluid particle with the consequence that microscopic 
molecular forces come into play and can no longer 
be neglected. In other words we are then out of the 
validity domain of the deterministic NSE and we 
would have to consider stochastic forces in addition 
to the deterministic external forces. It is therefore 
not yet settled, whether the phenomenon of hydro- 
dynamic turbulence is a manifestation of determin- 
istic chaos alone. 

As should be noted, the problem of finite time 
singularities cannot be overcome by some averaging 
recipe, because the existence of /(r) is connected to 
the uniqueness of solutions sufficient condition, 
and it may turn out to be also necessary. 

Similarly, in the case of compressible flows finite 
time singularities could not be excluded so far 0]. 
The proof or disproof of the existence of finite time 
singularities constitute one of the basic unsolved 
problems in the analysis of the NSE. In the inviscid 
case of the Euler equation, there is a general ar- 
gument for possible finite time singularities, see for 
instance Frisch |10|. From a direct numerical sim- 
ulation of the Euler equations, Grauer and Sideris 
pT[| recently reported on evidence for a singularity 
of the type as given in (37) with 7 = 1. 



Acknowledgement The author is indepted to 
M. Boudourides for making him aware of J. Serrin's 
contribution to the field. He is also thankful to A. 
Spille for a critical reading of the manuscript. 



5 



Appendix 

In the following it is sketched how the integral 
I(r) shows up in uniqueness proofs, see 0. At vari- 
ance with Q we do not adopt periodic boundary 
conditions. Let us assume there are two different 
solutions V and v' of the NSE (13). Then we define 
u := V — v' and obtain after subtracting the NSE 
for V and v' 



duk duk , duk , dvk d{p - p') 
— ^'"i^ = 

Ot OXi OXi OXi 
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+1/0 A Uk. 
(38) 

When this equation is scalarly multiplied by u and 
integrated over the volume 0, then, apart from the 
pressure term, the second and third terms of the left 
hand side vanish by the same argument used before 
in (31). With the abbreviation 



we can write 
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Now the viscosity term is estimated by means of the 
Poincare inequality p 



A 



1^1 > I HP 



(41) 



where / denotes the smallest distance between two 
parallel planes which just contain 0,. The B term 
is estimated by using the definition (36) and the 
Schwarz inequality as follows 
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dV \ukUi 



< 9||Dv||oo||uf . (42) 
One arrives at the ordinary differential inequality 
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which by Gronwall's lemma can be integrated to the 
final inequality 



|u(i)||2<||u(0)|pexp 



4uo 
' I 



2 t + 18I{t) 



(44) 



This result tells that, since the two supposed solu- 
tions V, v' possess the same initial conditions and 
therefore u(0) = 0, we have u(t) = for times 
t € (0, r) for which I{t) exists. This conclusion 
holds true also in the inviscid limit vq — > 0. 
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